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ABSTRACT 

The  combined  techniques  of  microbeam  x-ray  diffraction  and  Finite  Element  (FE) 
analysis  have  been  used  to  analyze  the  residual  stresses  resulting  from  an  interference 
operation  in  1080  steel.  Both  techniques  showed  the  resulting  stresses  to  be  tensile,  and 
accelerated  growth  rates  were  observed  in  fatigue  cracks  subsequently  grown  into  this 
residual  stress  field.  The  growth  rates  also  showed  a  clear  sensitivity  to  applied  stress  _ 
intensity  level.  Through-thickness  residual  stress  variations  as  predicted  by  FE  analysis 
help  explain  observed  crack  front  shapes. 

INTRODUCTION 

Whether  formed  during  processing  operations  or  accrued  in  service,  materials  will  often 
contain  residual  stresses.  In  structures  undergoing  fatigue  loading  these  residual  stresses 
can  alter  actual  loading  conditions  and  affect  the  growth  of  fatigue  cracks.  Thus  there 
have  been  a  number  of  smdies  undertaken  to  examine  residual  stress  effects  on  fatigue 
crack  growth  (see  [1]  for  a  review).  Although  compressive  residual  stresses  have  been 
almost  universally  found  to  decrease  growth  rates,  the  effect  of  tensile  residual  stresses  is 
not  as  clear.  Although  some  evidence  for  accelerated  growth  rates  in  the  presence  of 
tensile  stresses  exists  [2]  several  materials,  including  steels,  have  shown  these  stresses  to 
have  little  or  no  effect  on  stage  11  crack  growth.  The  lack  of  mean  (tensile)  stress 
sensitivity  in  the  Paris  growth  regime  has  been  linked  to  stress  fading  near  the  crack  tip, 
which  will  vary  depending  on  the  applied  stress  intensity  level  [3]. 

Comprehensive  analysis  of  residual  stress  effects  requires  stress  information  throughout 
the  body  of  the  a  material,  which  may  be  obtained  by  including  through-thickness  Finite 
Element  (FE)  stress  analysis  along  with  surface  sensitive  x-ray  diffraction  measurements. 
With  the  advent  of  greater  speed  and  flexibility  in  computer  systems  FE  modeling  of 
quantities  such  as  residual  stresses  and  plastic  deformation  is  taking  on  an  increased  role 
in  structural  analysis.  Confident  use  of  FE  to  model  real  problems,  however,  necessitates 
comparison  with  experimental  results  and  a  dearth  of  such  studies  exists.  Therefore  aFE 
code  was  used  to  model  the  development  of  tensile  residual  stresses  from  an  interference 
operation  and  meshed  to  compare  directly  near-surface  residual  stresses  with  measured 
x-ray  values.  Reasonable  agreement  was  obtained,  and  FE  predictions  of  through- 
thickness  residual  stress  variations  helped  explain  observed  crack  front  shapes. 

EXPERIMENTAL  PROCEDURE 

Specimens  and  FE  Mesh 

As  received  1080  steel  (12.4mm  thickness)  was  heat  treated  to  produce  a  fine  grained, 
spheroidized  microstructure.  The  resulting  average  ferrite  grain  size  and  cementite 

particle  size  was  5  and  O.Vum,  respectively,  with  a  13.1(6)%  cementite  volume  fraction 
as  determined  from  x-ray  powder  diffraction.  Both  sides  of  the  steel  were  electropolished 


to  remove  0.45mm  from  each  face,  resulting  in  smooth  surfaces  free  of  decarbeurized 
lavers  as  verified  by  hardness  measurements. 

romnart-tension  fCT)  specimens  were  machined  from  the  bulk  steel  as  per  ASp^I 
standards  [4],  and  residual  stresses  were  introduced  ahead  of  the  future  crack  planes  y 
driUin-  taperk  holes  and  then  expanding  the  holes  by  inserting  c^e-hardened  ^d 
tapered  stLl  pluss  at  0.95%  radial  interference  strain(fig  1(a)).  The  tapenng  was 
ZSs  Ji  to  Slow  uniform  through-thickness  pressure  to  be  applied  dunng  the 
interference,  but  led  to  an  -8%  larger  stress  on  the  bottom  surface  relative  to 
to  differences  in  plug  diameter.  Through-thickness  stress  vanations  ansing  from 
frictional  sliding  across  the  plug/plate  interface  were  also  expected,  u^ersconng 
importance  of  through-thickness  stress  information  provided  by  our  FE  model  m  additio 
to  surface-sensitive  x-rav  measurements.  'Control'  specimens  were  also  produce 
without  the  hole  and  with  the  hole  but  without  the  plug  for  direct  comparison  of  crack 

growth  rates. 

The  interference-fit  operation  was  modeled  along  the  future_(x-z)  crack  plane  by  FE  ^ 
shown  in  fig.  1(b).  Due  to  the  radial  symmetry  of  the  operation  this  x-z  slice  ®  ® 

axisymetrically,  with  the  center  of  the  plug  serving  as  the  axis  of  symmetry.  The 
interference  was  simulated  by  forcing  the  plug  to  expand  under  eigenstrains 
corresponding  to  the  applied  interference. 


(1) 

(2) 


e*=e-  =  5/(R  +raz)  =  0.03302/(3.5  +  0.0104z) 

=  o.oo3z 

where  5  is  the  constant  radial  misfit,  Ro  is  the  radius  of  the  top  of  the  plug  (at  z=0),  m  is 
the  taper  angle  in  radians  and  z  is  in  mm.  The  z-dependence  of  the  x  and  y  eigen.trams 
accounts  for  the  linear  variation  in  strain  (and  stress)  along  z  in  accordance  with_ 
experiment.  The  z-component  of  eigenstrain  simulates  the  fnctional  toce  expe^nced  by 
the^pluz/plate  interface.  Its  value  was  determined  empincally  by  matching  the  re¬ 
simulated  surface  displacements  with  expenmentd  values  me^ured  with  a  surface 
profilimeter  (not  shown).  The  FE  mesh  was  graded  towards  both  surfaces  and  the  ^ 

interface  -  with  the  smallest  surface  elements  5um  in  height  and  approximately  -OOiim  m 
width  near  the  interface.  This  allows  for  direct  companson  with  x-ray  data  taken  with  the 

200um  beam  as  the  1/e  penetration  depth  for  Cr  in  Fe  is  approximately  The 
resulting  mesh  was  comprised  of  3520  total  elements.  Matenal  parameters  for  1080  steel 

were  derermined  from  uniaxial  stress/strain  data  (ay=550MPa,  ^=2l5GPa)  ai^  mput 
into  the  Finite  Element  Analysis  Package  [5]  for  stress  analysis.  The  steel  s^jed  ^ 
aDoreciable  work  hardening  out  to  the  maximum  plastic  strain  ranges  computed  by  FE 
aS  so  was  modeled  as  elastic-perfectly  plastic,  while  the  much  harder  plug  was  taken  t 
be  elastic  throughout  with  the  same  modulus  as  the  plug. 

Fatigue  Testing 

The  CT  specimens  were  notched  with  a  jeweler’s  saw  and  pre-cracked  by  S  ^ 

tensile  load,  which  was  then  decreased  in  steps  to  grow  the  cracks  into  the  stable  sta„e  U 

regime  at  the  desired  AK  values.  These  stage  E  fatigue  cracks  were  cycled  at  R=0.1  and 
co^'nstant  AK  at  levels  ranging  from  19.5-25  MPaV^.  Through-thickness  crack  ^owth 
was  monitored  by  measuring  specimen  c°rriphance  with  a  clip  gage  and^^^^  Y 
checked  with  an  optical  microscope  until  sample  failure  into  the  plu^  int 


X-ray  Microbeam  Stress  Analysis 

Diffraction  methods  are  widely  used  for  stress  analysis  due  to  their  relatively  non¬ 
destructive  nature  and  ability  to  measure  the  entire  stress  tensor.[6]  This  analysis  is  based 
on  measuring  average  lattice  strains  from  diffraction  peak  shifts,  which  are  converted  to 
stresses  using  measured  elastic  constants.  In  multiphase  materials  the  stresses  can  be 
measured  in  each  phase  allowing  separation  into  macrostress  and  microstress 
components. 

Here  we  used  triaxial  stress  analysis  to  measure  both  in-  and  out-of-plane  stress 
components  via  a  least-squares  procedure[7],  and  separated  stresses  into  macro-  and 
microstresses  in  each  phase  using  the  measured  cementite  volume  fraction  and  equations 
of  microstress  equilibrium  [8]. 

In  order  to  measure  stress  gradients  generated  near  the  plug  interface  an  x-ray  microbeam 

of  2  lOum  diameter  was  used.  Since  beam  intensity  is  proportional  to  area  such  a 
microbeam  is  relatively  weak  -  thus  we  designed  our  system  to  maximize  flux  and 
efficiency.  A  tapered  glass  capillary  was  developed  to  focus  x-rays  produced  by  a  12 
kW  rotating  anode  source  with  a  Cr  target.  The  capillary  provides  a  threefold  increase  in 
intensity,  and  comparable  increase  in  divergence,  as  compared  to  an  equivalent  pinhole 
collimator.  The  Fe  2 1 1  peaks  were  measured  with  a  position  sensitive  detector  while  a 
solid-state  detector  was  used  for  the  weaker  Fe3C  250  peaks  in  order  to  increase  the 
peak/background  ratio.  Further  details  of  the  experimental  setup  and  measurements  are 
given  in  [9]. 

RESULTS  AND  DISCUSSION 

X-ray  and  FE  Residual  Stresses  from  Interference  Fit 

Comparisons  between  the  measured  and  simulated  residual  stresses  formed  in  the  near- 

surface  regions  (<5jim)  of  the  plate  are  shown  in  fig2(a,c).  The  deviatoric  stresses  'Cyy 
are  shown  since  these  will  most  strongly  affect  the  crack  during  mode  I  loading.  Plastic 
deformation  as  evidenced  from  diffraction  line  broadening  was  measured  out  to  4mm 
from  the  interface.  Within  this  plastic  region  the  cementite  takes  up  a  higher  fraction  of 
the  load  than  the  ferrite,  leading  to  microstresses  in  addition  to  a  tensile  macrostress.  The 
FE  method  qualitatively  reproduces  these  measured  macrostresses  on  both  surfaces  of  the 
plate,  though  it  overestimates  their  magnitude  near  the  interface. 

By  qualitatively  matching  the  measured  surface  macrostresses  FE  predictions  of 

throughout  the  crack  plane,  fig.2(c),  can  be  used  with  greater  confidence  to  predict  z- 
variations  in  residual  stress  and  therefore  crack  front  shape.  Two  major  trends  seen  from 
the  FE  predictions  are  that  the  tensile  stresses  are  greatest  at  the  surfaces  of  the  plate  and 
that  the  frictional  sliding  simulated  along  the  plug/plate  interface  has  led  to  especially  high 
tensile  stresses  near  the  top  of  the  plate  relative  to  the  bottom  portion.  Thus  as  the  crack 
approaches  the  interface  growth  is  expected  to  be  accelerated  at  the  surface  relative  to  the 
bulk,  especially  towards  the  top  surface  of  the  plate. 

Fatigue  Crack  Growth  Rates  and  Crack  Front  Shape 

Fatigue  cracks  grown  into  the  tensile  stress  fields  were  significantly  affected  by  these 
stresses  as  shown  in  fig  3(a).  In  this  figure  the  propagation  rates  at  each  AK  were 

normalized  by  the  rates  of  control  samples  grown  at  equal  AK,  thus  limiting  any 
differences  in  growth  rates  to  residual  stress  effects.  The  tensile  residual  stresses  increase 


crack  growth  rates  by  up  to  two  orders  of  magnitude  near  the  plug.  Furthermore,  the 
distance  at  which  the  crack  begins  accelerating  shows  a  strong  dependence  on  AK.  At 
AK=  1 9.5  MPaA/m  the  crack  begins  to  accelerate  at  a/c=0.74  (6.5mm  from  the  plug/plate 

interface)  while  as  AK  is  increased  the  acceleration  is  progressively  retarded  to  distances 
closer  to  the  interface.  This  implies  that  the  crack  becomes  more  'sensitive'  to  residual 

stresses  at  lower  applied  AK.  Once  the  acceleration  has  begun  at  a  given  AK,  however,  it 

proceeds  at  a  higher  rate  with  increasing  AK,  which  may  be  associated  with  a  transition 
from  stage  n  to  catastrophic  stage  m  growth  under  the  combination  of  high  applied  and 
residual  tensile  stresses.  To  ensure  that  these  effects  are  not  geometrically  based  we  grew 
a  crack  into  a  sample  containing  a  hole  without  the  plug  and  found  crack  growth  to  be 
unaffected  (insert).  This  shows  that  the  residual  stresses,  father  than  geometrical  effects, 
are  causing  the  observed  acceleration. 


Fractographic  analysis  on  the  sample  grown  at  AK=22.5  MPaV m  into  the  plug  stress 
field,  fig.3(b),  illustrates  the  effect  of  the  through-thickness  stress  variations  as  predicted 
by  FE.  While  crack  fronts  in  control  samples  showed  bowing  towards  the  interior  due  to 
the  plane  strain  (high  driving  force)  conditions  relative  to  the  plane  stress  conditions  at  the 
surface,  the  crack  fronts  in  the  plug  sample  are  quite  different.  At  3mm  from  the  interface 
(position  A  in  the  figure)  the  crack  front  is  rather  linear,  evidence  that  the  higher  surface 
stresses  (relative  to  the  bulk)  have  accelerated  surface  growth,  while  closer  to  the  plug 

(position  B)  the  crack  has  grown  asymmetrically  towards  the  position  of  maximum  cr^ 
as  predicted  by  FE  (fig.2(b)).  ^ 

CONCLUSIONS 

1.  An  interference  fit  in  1080  steel  resulted  in  a  complex  residual  stress  field,  primarily 
consisting  of  tensile  macrostresses,  which  was  characterized  by  microbeam  x-rav 
diffraction  and  comparative  Finite  Element  (FE)  analysis. 

2.  These  tensile  stresses  significantly  accelerated  fatigue  crack  growth  rates,  with 

increasing  crack  sensitivity  to  the  residual  stresses  observed  with  decreasin'^  applied 
stress  intensity.  °  ^ 

3.  Through-thickness  residual  stress  variations  as  predicted  by  FE  analysis  help  explain 
observed  crack  front  shapes  and  thus  highlight  the  usefulness  of  FE  in  conjunction  with 
comparative  expenmental  analysis  to  understand  residual  stress  effects  in  fati-nie 
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Figure  1:  Geometry  for  (a)  CT  specimen  and  (b)  corresponding 
Finite  Element  slice,  with  element  grading  and  displacement  boundary 
conditions  shown  on  mesh  edges.  Unless  given,  all  dimensions  in  mm. 
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